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In  the  course  of  our  previous  studies  of  elastic  waves  generated  by 
nuclear  explosions,  a  number  of  questions  have  arisen  which  are  difficult 
to  answer  through  the  analysis  of  data  from  individual  seismic  stations. 
The  direction  of  approach  of  various  phases,  their  apparent  velocity,  and 
the  degree  of  their  spatial  coherence  are  a’’  matters  which  are  not 
completely  rede-stoou  for  data  recorded  at  small  distances  from  explosions, 
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This  lack  of  understanding  hinders  the  interpretation  of  the  data  and  the 
study  of  the  explosive  source  mechanism.  Such  considerations  have  led  us 
to  use  a  small  array  of  seismic  stations  to  record  the  elastic  wave  field 
from  a  number  of  explosions. 

A  seismic  array  consisting  of  9  stations  arranged  as  three  nested 
triangles  with  minimum  spacing  of  100  m  and  maximum  spacing  of  400  m 
was  used  to  record  the  two  NTS  explosions,  LIPTAUER  in  Yucca  Valley  on 
3  April  1980  and  COLWICK  on  Pahute  Mesa  on  26  April  I98O.  The  distance 
of  the  array  from  the  explosion  was  1.9  km  for  LIPTAUER  and  6.0  km  for 
COLWICK. 

So  far  our  analysis  has  concentrated  on  the  COLWICK  data.  Preliminary 
results  show  that  most  of  the  energy  in  the  first  5  sec  of  the  seismograms 
is  propagating  at  an  azimuth  and  velocity  which  are  consistent  with 
generation  at  or  very  the  explosive  source.  The  spatial  coherence 

decays  by  a  factor  ofie"  at  distances  that  range  between  AOO  and  1200  m, 
depending  upon  direction  and  component  of  motion.  In  a  direction  parallel 
to  the  wavefront,  the  vertical  and  radial  components  are  much  more  coherent 
than  the  transverse  component,  while  all  of  the  components  have  comparable 
coherence  in  a  direction  perpendicular  to  the  wavefront.  Coherence  generally 
decreases  as  frequency  increases. 

In  addition  to  aiding  our  interpretation  of  elastic  waves  generated 
by  an  explosion,  these  array  data  should  provide  useful  constraints  on  the 
possible  mechanisms  by  which  elastic  waves  can  be  scattered  by  inhomogeneities 
near  the  surface  of  the  earth.  ' 

A  large  body  of  surface  waVe  data  recorded  at  regional  distances  in  the 
Cal ifornia-Nevada  region  is  being  analysed  to  obtain  estimates  of  velocity 
and  attenuation  in  the  crust.  The  observed  surface  waves  require  a 
considerable  amount  of  processing  in  order  to  extract  the  necessary  information 
Both  the  phase-matched  filter  and  moving-window  analysis  schemes  have  been 
used  with  the  former  proving  to  be  superior.  It  appears  that  excellent 
dispersion  and  amplitude  data  can  be  extracted  from  the  seismograms  in  the 
period  range  5  to  50  sec. 

For  the  sake  of  completeness,  a  report  is  included  which  was  prepared 
in  response  to  a  specific  DARPA  request.  This  report  is  a  review  of 
"Near-source  effects  on  P  waves". 

A  catalog  of  digital  data  has  been  assembled  whici,  cc'^ains  much  of 
the  data  collected  and  analysed  in  DARPA  supported  resec  of  previous 
years.  A  preliminary  version  of  the  catalog  index  for  the  period  1965-1971 
is  included  in  this  report. 
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than  the  transverse  component,  while  all  of  the  components  have  comparable 
coherence  in  a  direction  perpendicular  to  the  wavefront.  Coherence  generally 
decreases  as  frequency  increases. 

In  addition  to  aiding  our  interpretation  of  elastic  waves  generated 
by  an  explosion,  these  array  data  should  provide  useful  constraints  on  the 
possible  mechanisms  by  which  elastic  waves  can  be  scattered  by  inhomogeneities 
near  the  surface  of  the  earth. 

A  large  body  of  surface  wave  data  recorded  at  regional  distances  in  the 
Ca 1 i forn i a-Nevada  region  is  being  analysed  to  obtain  estimates  of  velocity 
and  attenuation  in  the  crust.  The  observed  surface  waves  require  a 
considerable  amount  of  processing  in  order  to  extract  the  necessary  information. 
Both  the  phase-matched  filter  and  moving-window  analysis  schemes  have  been 
used  with  the  former  proving  to  be  superior.  It  appears  that  excellent 
dispersion  and  amplitude  data  can  be  extracted  from  the  seismograms  in  the 
period  range  5  to  50  sec. 

For  the  sake  of  completness,  a  report  is  included  which  was  prepared 
in  response  to  a  specific  DARPA  request.  This  report  is  a  review  of 
"Near-source  effects  on  P  waves". 

A  catalog  of  digital  data  has  been  assembled  which  contains  much  of 
the  data  collected  and  analysed  in  DARPA  supported  research  of  previous 
years.  A  prel imi nary  version  of  the  catalog  index  for  the  period  1965-1971 
is  included  in  this  report. 


T 


I 


1 


1  .  Summary 

In  the  course  of  our  previous  studies  of  elastic  waves  generated  by 
nuclear  explosions,  a  number  of  questions  have  arisen  which  are  difficult 
to  answer  through  the  analysis  of  data  from  individual  seismic  stations. 

The  direction  of  approach  of  various  pha:>es,  their  apparent  velocity,  and 
the  degree  of  their  spatial  coherence  are  all  matters  which  are  not 
completely  understood  for  data  recorded  at  small  distances  from  explosions. 
This  lack  of  understanding  hinders  the  interpretation  of  the  data  and  the 
study  of  the  explosive  source  mechanism.  Such  considerations  have  led  us 
to  use  a  small  array  of  seismic  stations  to  record  the  elastic  wave  field 
from  a  number  of  explosions. 

A  seismic  array  consisting  of  9  stations  arranged  as  three  nested 
triangles  with  minimum  spacing  of  100  m  and  maximum  spacing  of  AOO  m 
was  used  to  record  the  two  NTS  explosions,  LIPTAUER  in  Yucca  Valley  on 
3  April  1980  and  COLWICK  on  Pahute  Mesa  on  26  April  I98O.  The  distance 
of  the  array  from  the  explosion  was  1.9  km  for  LIPTAUER  and  6.0  km  for 
COLWICK. 

So  far  our  analysis  has  concentrated  on  the  COLWICK  data.  Preliminary 
results  show  that  most  of  the  energy  in  the  first  5  sec  of  the  seismograms 
is  propagating  at  an  azimuth  and  velocity  which  are  consistent  with 
generation  at  or  very  near  the  explosive  source.  The  spatial  coherence 
decays  by  a  factor  of  e  ^  at  distances  that  range  between  400  and  1200  m, 
depending  upon  direction  and  component  of  motion.  In  a  direction  parallel 
to  the  wavefront,  the  vertical  and  radial  components  are  much  more  coherent 
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I  I  .  Near  field  array  studies 


no  tivaticn 


From  Pahute  Mesa  NTS  explosions  anomalous  features  in  the  near 
field  accelerograms  have  been  noticed  ^in  previous  field  experiments. 
Significant  transverse  accelerations  are  observed  coincident  with  the 
initial  P  wave  arrival.  It  was  not  clear  whether  this  transverse  faction 
was  coherent  for  any  considerable  areal  extent, or  was  confined  to  local¬ 
ized  sites.  Our  previous  experiments  with  widely  spaced  stations  were 
not  suitable  to  answer  this  question.  The  energy  of  the  transverse 
seismogram  is  comparable  to  the  vertical  seismogram  and  transverse  peak 
accelerations  often  exceed  the  vertical  or  radial  peak  accelerations. 
Due  to  the  sensitivty  of  the  transverse  components  to  a  deviatoric 
(non-explosive)  source  component,  it  was  desirable  to  understand  the 
nature  of  the  waves  propagating  with  an  apparent  transverse  component, 
and  to  test  the  hypothesis  that  large  transverse  motions  were  due  to 
off-azimuth  SV  waves  as  opposed  to  on-azimuth  SH  waves. 

Design 

A  small  array  of  closely  spaced, broadband ,  large  dynamic  range, 
three-component  accelerometers  was  considered  to  address  the  questions 
on  spatial  coherence  and  wave  properties.  As  with  previous  experiments, 
digital  event  recorders  with  three-component  packages  of  force-balanced 
servo-accelerometers  were  used.  Ideally,  we  needed  to  measure,  1. 
direction  of  arrival,  2.  apparent  velocity  and,  spatial  coherence. 
The  array  had  to  be  easily  and  quickly  deployed,  to  be  fault  tolerant, 
to  be  limited  to  nine  stations,  to  afford  common  timing  and  to  preclude 
spatial  aliasing  while  providing  sufficient  resolution  in  two  dimen¬ 
sional  slowness  space.  A  nested  set  of  three  triangles  was  selected. 

Deployment  was  made  for  two  shots,  Liptauer  and  Colwick  at  1.9  and 
6.0  kilometers  respectively.  Liptauer  was  an  equivalent  local  magnitude 
4.7  explosion  located  in  Yucca  Valley.  Colwick  was  an  equivalent  local 
magnitude  5.5  explosion  on  Pahute  Mesa.  At  Liptauer  a  common  triggering 
of  the  array  was  attempted  to  provide  relative  timing.  At  Colwick  all 
event  recorders  were  allowed  to  trigger  individually  and  all  records 
were  time  shifted  assuming  a  given  P  wave  velocity.  At  each  site  one 
out  of  nine  recorders  malfunctioned,  leaving  eight  stations  arranged  as 
in  Figures  2  and  3-  The  wave  number  impulse  responses  are  shown  in  Fig¬ 
ures  4, and  5.  Resolution,  for  each  array,  was  about  1  cycle/km  at  half 
maximum  response  and  the  spatial  Nyquist  wavenumber  was  about  6  cycles/ 
km  in  all  directions.  It  was  possible  to  make  crude  slowness  spectral 
estimates  in  the  1  to  4  Hertz  band.  The  spatial  decay  of  coherence 
could  be  observed  in  the  longitudinal  as  well  as  the  transverse 


directions.  For  any  frequency  band  the  rate  at  which  coherence  between 
stations  decays  with  increasing  separation  could  be  measured  in  two 
dimension. 

Deployment  of  the  array  required  four  to  six  man  days  to  survey, 
implace,  and  test.  Removal  of  the  equipment  required  between  one  and 
two  man  days,  less  travelling.  The  arrays  were  deployed  with  in  a  week 
prior  to  the  event  and  external  power  was  provided  by  automobile  twelve 
volt  batteries  v/ith  sufficient  amp-hours  to  power  the  accelerometers  and 
digital  event  recorders  for  one  week.  Common  triggering  required  laying 
a  tv/isted  pair  of  wires  between  stations  to  provide  a  common  slave  sig¬ 
nal.  The  trigger  signal  logic  required  three  simultaneous  trigger  sig¬ 
nals  from  three  of  the  nine  event  recorders.  Aproximately  one  half 
second  of  pre-event  memory  was  availible  at  200  samples  per  second  for 
all  three  componeiits.  Records  are  reJughly  fifty  seconds  long. 
Antialias  filters  v/ere  either  at  25  or  50  hertz. 

Liptauer 

The  traces  from  the  Liptauer  array  are  seen  in  Figures  6,  7, and  8. 
The  ratio  of  the  vertical  to  radial  amplitudes  for  both  P  and  S  waves 
and  the  short  duration  of  the  seismograms  suggest  a  half  space  model  may 
fit  the  data  well.  The  largest  peak  transverse  signal  rivals  the  radial 
component  and  often  exceeds  it  at  individual  stations.  The  large 
transverse  pulse  at  1.7  sec  is  visually  coherent.  The  slight  emergent 
transverse  motion  coincident  with  the  P  wave  is  coherent  as  well. 

Colwick 

The  traces  from  the  Colwick  array  are  seen  in  Figures  9,  10, and  11. 
The  largest  peak  acceleration  was  18+/-3  ^  g  on  the  radial  component. 
The  large  vertical-to-radial  P  wave  ratio  and  radial-to-vertical  S  wave 
ratio  suggest  velocity  gradients  to  increase  the  angles  of  incidence. 
The  lengths  and  complexity  of  the  seismograms  also  indicate  that  a  more 
complicated  structure  may  be  necessary  for  Pahute  Mesa  than  is  needed 
for  Yucca  Valley.  Visually,  the  Colwick  seismograms  show  less  coher¬ 
ence,  possibly  due  to  the  presence  of  more  high  frequencies. 

Processing 

Programs  were  developed  to  perform  most  of  the  desired  processing 
upon  a  captive  minicomputer.  The  remainder  of  the  computing  was  per¬ 
formed  on  a  CDC  6400.  The  individual  event  recorder  tapes  are  played 
into  the  minicomputer,  the  data  are  reformated  and  stored  on  disk  and 
tape.  Processing  includes  time  domain  stacking,  spectra,  high  resolu¬ 
tion  and  convential  wave-number  analysis,  and  cross-correlations.  At 
present  processing  continues  upon  both  data  sets  while  the  analysis  for 
the  Colwick  array  is  more  complete.  Discussion  will  stress  the  results 
obtained  to  date  from  the  Colwick  array. 

1  .  Stacking 

From  the  wave  n'omber  impulse  responce  we  can  estimate  that  the 
resolution  of  the  Colwick  array  at  5  Hz  is  +/-  0.2  sec/kra.  Stacking  of 
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the  seisinogram;5  was  done  for  sloivnesses  from  1.2  to  -1.2  sec/km  at  steps 
of  0.2  sec/km  on  azimuth  with  the  shot.  The  results  indicate  that  the 
majority  of  the  coherent  energy  propogates  from  the  shot  at  slownesses 
less  than  1  sec/km  and  that  the  initial  body  waves  all  arrive  with 
slownesses  less  than  0.25  sec/km,  with  respect  to  the  initial  P  wave. 
Stacking  of  the  seismograms  can  be  seen  tu  deemphas’ the  high  frequen¬ 
cies.  Unless  these  high  frequencies  propagate  at  significantly  higher 
apparrent  velocities,  there  must  be  a  loss  of  coherency  with  increasing 
frequency.  In  an  attempt  to  measure  this  phenomenon,  filtered  two- 
station  cross-correlations  were  calculated.  The  results  are  presented 
in  section  3.  The  large  transverse  pulse  at  1.5  sec  stacks  up  0.2 
sec/km  slower  than  the  P  wave.  Some  later  arrivals  have  been  found  to 
stack  up  for  off  azimuth  slownesses.  A  stategy  that  has  been  useful  for 
detecting  scattered  waves  is  to  stack  suspicious  off-azimuth  peaks  that 
occur  in  the  wave  number  spectra.  , 

2.  Wavenumber  spectra 

Windows  were  selected  for  the  calculation  of  cross-spectra  to  com¬ 
pute  slowness  spectra  at  selected  frequencies.  Both  conventional  and 
"high  resolution"  methods  were  used.  The  beam  forming  or  conventional 
method  retains  the  dynamic  range  of  the  data  but  is  cluttered  by  side 
lobe  excitation.  The  partial  removal  of  side  lobes  by  the  high  resolu¬ 
tion  or  maximum  liklihood  method  is  convienant  for  these  arrays  (  Capon, 
Greenfield,  and  Kolker,  1967;  Capon,  1969;  and  Capon,  1970  ). 

The  preliminary  analysis  of  the  Colwick  and  Liptauer  array  data 
show  that  coherent  energy  from  the  source  azimuth  is  dominant  on  all 
three  components  for  the  first  5  seconds.  In  Figures  13>  and  14  we  see 
the  Colwick  array  high  resolution  slowness  spectra  at  3*7  Hz  for  the 
first  5  seconds  of  data.  The  origin  of  the  plots  corresponds  to  the 
slowness  of  the  initial  P  wave  arrival.  The  convention  is  for  the  wave 
number  vector  to  project  in  the  direction  from  which  the  arrival  comes. 
The  shot  was  in  the  direction  of  the  negative  Y  axis.  The  vertical  com¬ 
ponent  acceleration  f Figure  13)  shows  a  prominent  peak  0.25  sec/km  fas¬ 
ter  than  the  Initial  ?  wave  and  hence  in  the  direction  of  the  plus  Y 
axis.  This  confirms  the  impression  derived  from  time  domain  stacking 
that  some  body  waves  are  very  fast.  The  remainder  of  the  energy  in  the 
vertical  component  is  distributed  at  slownesses  less  than  0.7  sec/km 
with  respect  to  the  P  wave  on  the  azimuth  of  the  source. 

The  large  transverse  pulses  are  not  off-azimuth  SV  waves  to  within 
the  resolution  of  the  arrays.  The  large  transverse  pulse  at  1.5  sec 
observed  at  the  Colwick  array  coincides  with  a  well  defined  peak  in 
slowness  space  0.2  sec/km  slower  than  the  P  wave,  as  seen  in  Figure  14. 
This  is  consistent  with  propagation  as  a  direct  S  wave  from  the  source. 
The  remainder  of  the  coherent  energy  in  the  transverse  component  is  dis¬ 
tributed  between  0.0  and  1.0  sec/km,  with  respect  to  the  initial  P  wave, 
and  on  azimuth  with  the  shot.  During  these  first  5  seconds,  the 
transverse  component  shows  no  indication  of  contamination  by  coherent 
arrivals  faster  than  the  initial  P  wave.  Windows  later  in  the  records 
have  revealed  significantly  off-  and  back-azimuth  arrivals  which  stack 
up  in  the  time  domain  as  distinct  arrivals.  However,  these  scattered 
packets  of  energy  are  not  observed  when  the  first  five  seconds  of  data 
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are  analyzed  as  a  block. 

3.  Two-Station  Cross-Correlations 


Both  time  domain  stacking  and  wavenumber  analysis  utilize  tne  whole 
array  and  emphasize  the  coherent  parts  of  the  seismogram.  In  order  to 
explore  the  incoherent  part  of  the  seismogram  vie  formed  two  station 
cross-correlations.  Out  of  eight  stations  we  can  make  28  pairs  of  sta¬ 
tions.  The  Colwick  array  was  arranged  to  order  data  into  transversely 
and  longitudinally  separated  station  pairs.  This  provides  a  comparison 
of  repeatable  ground  motion  across  a  wave  front  as  well  as  along  the 
direction  of  propagation. 

• 

The  cross-correlation  function  serves  as  an  unbiased  estimate  of 
the  spatial  covariance  (Weiner,  1949;  Bretherton  and  McWilliams  1980). 
Furthermore,  filtered  cross-correlations  yield  information  about  the 
frequency  dependence  of  coherence.  Analysis  of  the  spatial  decay  of 
coherence  can  serve  to  test  theories  of  wave  propagation  in  randomly 
heterogeneous  media.  Such  formalisms  for  scalar  wave  propagation 
predict  differences  between  longitudinal  and  transverse  spatial  covari¬ 
ance  functions  (Chernov  1 960, Ishimaru  1979).  Analysis  has  progressed  in 
this  direction  for  the  Colwick  array  data.  We  have  concentrated  the 
processing  upon  the  first  5  sec  block  of  the  data  for  a  test  of  metho¬ 
dology. 


The  two-station  cross-correlations  for  all  three  components  were 
considered  to  be  functions  of  station  separation  and  separation  direc¬ 
tion.  For  stations  i  and  j  at  locations  r^  and  r.  the  cross-correlation 
.shown  in  Figure  15,  was  modeled  as 


In  particular,  a  fair  fit  to  the  data  was  obtained  by  assuming  a  sim¬ 
ple  exponential  decay. 


R.  .  =  exp(-a  Ix.-x.I-a  L  . -y  .  j ) 

W'ith  the  X  axis  aligned  in  the  transverse  direction,  and  the  y  axis  in 
the  longitudinal  direction.  Because  the  array  had  only  13  independant 
pairs  of  1 x. -x  .  j , 1 y . -y  .  j  out  of  28  pairs  of  stations,  consistency  of  the 
data  could  fee  examined^  In  Figure  16  we  see  the  values  of  the  data  con¬ 
toured  by  the  fit  to  a  ,a  for  each  component.  The  RMS  residuals  for  the 
three  components  indicate^a  fair  fit. 


1  .0/  a 

X 

1  .0/ 

rms 

vertical 

1465  meters 

663  meters 

.04 

rad ial 

1960 

471 

.07 

transverse 

429 

477 

.11 

residual 


For  comparison,  the  two  station  cross-correlations  are  plotted  in 
Figure  17  versus  the  station  separation.  This  one-dimensional 
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representation  shov/e  much  scatter.  We  see  throe  curves  fo”  each  com¬ 
ponent.  Tliese  represent  cross  sections  through  the  t«ro-dimrnsionnl  sur¬ 
face  tha'’  was  fit  to  the  data.  For  a  perfect  fit  to  the  data  these 
curves  would  span  the  the  individual  "easurements  plotted  in  Figure  17. 
The  data  suggest  we  can  rank  the  components  in  increasing  sensitivity  to 
heterogenei ty ,  for  an  explosion  source,  as  vertical,  radial,  and 
t  ra ns verse . 

The  loss  of  coherence  along  the  wave  front  (transverse  separara- 
tion)  for  the  transverse  component  is  of  the  same  order  as  in  the  direc¬ 
tion  of  propagation.  T  o  radial  and  vertical  components  show  more 
coherence  across  the  the  wave  front  than  in  the  dirction  of  propagation 
at  the  same  separation.  It  must  be  noted  that  the  loss  of  coherence  in 
the  longitudinal  direction  due  to  a  changing  deterministic  seismogram 
with  distance  has  not  been  removed.  It  mifcl  also  be  emphasized  that 
this  represents  the  whole  of  the  first  5  seconds  of  data,  including  both 
coherent  signal  seen  in  the  wave  number  spectra  and  incoherent  energy 
that  does  not  stack  in  beam  forming.  This  incoherent  energy  can  be  con¬ 
sidered  a  form  of  signal  generated  noise  produced  by  an  interaction  of 
the  geologic  heterogeneity  with  the  coherent  waves  propagating  througr' 
an  average  homogeneous  structure. 

Additional  information  exists  in  the  frequency  dependence  of  cuher- 
ence  (Knopoff  and  Hudson,  1964  and  1967;  Hudson  and  Knopoff,  1966;  and 
Karal  and  Keller,  1964).  Each  component  of  ground  motion  contains  scat¬ 
tered  energy  contributed  by  all  raodes  of  propogation.  A  source  with 
small  deviatoric  part  produces  dominating  radial  and  vertical  com¬ 
ponents.  The  data  suggest  the  incoherent  portion  of  the  transverse  com¬ 
ponent  of  ground  motion  is  dominated  by  scattc^red  energy  from  the  P  and 
SV  modes  from  the  explosive  source.  The  remaining  transverse  coherent 
energy  is  propagating  SH  waves  from  a  small  deviatoric  source.  Due  to 
the  larger  initial  excitation  of  the  P  and  SV  nodes  they  contain  ,on  the 
average,  proportionately  less  incoherent  scattered  waves.  It  has  been 
shown  (Knopoff  and  Hudson  1967  )  that  the  P-to-S  and  S-to-S  conversions 
are  dominant  over  ?-to-P  and  S-to-P  conversions.  The  P  and  SV  modes 
of  propagation  will  supply  incoherent  enei’gy  to  the  SH  modes  of  propoga- 
ticn  in  much  larger  proportion  than  will  the  SH  modes  provide  scattered 
energy  to  the  P  and  SV  modes  of  propogation.  The  resultant  tendency 
will  be  for  the  incoherent  fraction  of  the  transverse  component  to  be 
larger  for  the  radial  and  vertical  fractions  of  incoherent- to-coherent 
energy.  Frequency  dependence  of  the  coherence  for  all  three  components 
will  be  helpful  to  estimate  inhomogeneity  in  the  medium. 

An  estimation  of  the  frequency  dependence  of  coherence  was 
attempted  by  filtering  the  cross-correlation  functions.  Seven  pairs  of 
stations  were  selected  that  gave  either  purely  transverse  separation  or 
purely  longitudinal  separation.  In  order  to  remove  the  frequency  depen¬ 
dence  of  the  source,  the  cross-spectra  were  normalized.  The  cross¬ 
spectra  were  then  low-pass  filtered  at  2.5,  5*0,  10.,  and  20.  hertz. 
The  inverse  Fourier  transform  was  applied  yielding  a  filtered  cross- 
correlation.  We  see  an  example  for  stations  1  and  4  in  Figure  18  .  We 
can  interpret  the  decreasing  correlation  with  in  increasing  bandwidth  as 
showing  91  ^  of  the  bandwidth  below  5  hertz  is  coherent  while  only  64? 
coherent  below  20.  hertz. 
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^or  the  three  components  the  maximum  correlations  are  plotted  in 
Fi^^ure  19  versus  the  corner  fi-enuency  of  the  filter.  The  station  pairs 
are  separated  into  two  populations  purely  transverse  or  purely  longitu¬ 
dinal  reparations.  All  three  components  show  a  sharp  fall-off  with 
increasing  corner  frequency.  Although  the  vertical  and  radial  com¬ 
ponents  show  the  two  populations  as  separated,  the  transverse  comnonent 
can  not  distinguish  between  the  two  populations.  Further  analysis  will 
include  bandpass  filtering  of  the  cross-correlations  as  well  as  detailed 
analysis  of  selected  shorter  windows.  The  exact  nature  of  the  energy  on 
the  transverse  component  prior  to  the  arrival  of  a  causal  3  wave  will 
require  such  a  window. 

Conclusions 

The  analysis  of  the  Colwick  array  datS  is  still  incomplete,  and 
much  remains  to  be  done  with  the  Liptauer  experiment  data.  Comparison 
of  the  two  sites  will  be  useful  for  a  test  of  the  methodology  as  well  as 
a  test  of  the  two  different  geologic  settings.  Preliminary  analysis  fo'’ 
both  arrays  has  demonstrated  that  a  significant  fraction  of  the 
transverse  energy  consists  of  propagating  SH  waves  arriving  at  apparent 
velocities  and  tines  consistant  with  a  deviatoric  (non-explosive)  com¬ 
ponent  in  the  source  region. 

The  arrays  of  closely  spaced  stations  provide  a  means  to  investi¬ 
gate  the  spatial  coherence  of  waves  in  a  complex  geologic  setting.  Our 
confidence  in  the  use  of  accelerograms  to  study  the  source  has  been 
increased.  The  persistence  of  the  coherent  energy  of  the  early  waves 
arriving  at  the  expected  azimuth  supports  the  use  of  laterally  homogene¬ 
ous  Green's  functiDns  in  source  studies. 

At  the  same  tine, the  array  data  provide  definition  of  the  statisti¬ 
cal  coda  contained  In  the  seismogram.  The  sensitivity  of  the  components 
of  ground  notion  to  lateral  heterogeneity  for  those  explosion  sources 
has  been  measured.  The  radial  and  vertical  components  were  found  to 
have  comparable  amounts  of  incoherent  energy  increasing  proportionately 
with  frequency.  Tne  conerence  along  the  wave  front  for  these  components 
was  found  to  be  excellent  in  a  tondwidth  up  to  10  Hz  at  ranges  exceeding 
the  size  of  the  arrays.  The  transverse  component  of  the  ground 
acceleration  was  found  to  differ  from  the  other  two  components  in  the 
proportion  of  scattered  energy.  In  terms  of  equivalent  separation  for 
given  loss  in  waveform  coherence,  a  mislocation  of  AOO  meters  at  6  km  in 
a  transverse  direction  would  affect  the  SH  seismogram,  by  the  same 
amount  as  would  a  similar  positional  error  in  the  radial  distance  from 
the  event.  For  radial  or  vertical  seismograms,  a  mislocation  in  the 
racial  direction  of  400  meters  would  affect  the  seismogram  the  same 
amount  as  a  positiond  error  in  the  azimuthal  direction  of  as  much  as 
1.5  km.  We  may  think  ef  tHpfs  incoherent  energy  as  signal  generated 
noise  and  treat  it  accordingly  for  the  purposes  of  inversion  procedures. 

The  remaining  work  to  be  done  will  req.ire  an  estimate  of  the 
coherent- to-incoherent  signal  ratios  for  the  various  components  in  the 
appropriate  geologic  settings.  This  is  in  a  very  broad  sense  the  deter- 
nxr.^i-ion  of  the  statistics  for  the  stochastic  Green's  function  ( Adomain 
1964). 
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Figure  Captions 

Figure  1.  Locations  of  the  NTS  explosions  which  were  studied  and  the 
accelerometer  arrays. 

Figure  2.  Arrangement  of  the  accelerometer  array  for  the  explosion  LIPTAUER. 

Figure  3-  Arrangement  of  the  accelerometer  array  for  the  explosion  COLWICK. 

Figure  Wavenumber  impulse  response  for  the  LIPTAUER  array. 

Figure  5.  Wavenumber  impulse  response  for  the  COLWICK  array. 

Figure  6,  7,  8.  13  sec  of  individual  traces  from  the  LIPTAUER  array. 

Station  numbers  are  on  the  left.  The  vertical  bars  are  0.1  g. 

Figures  9>  10,  11.  5.5  sec  of  individual  traces  from  the  COLWICK  array. 

Station  numbers  are  on  the  left.  The  vertical  bars  are  0.1  g. 

Figure  12.  Averaged  acceleration  spectra  from  the  first  5  sec  at  the 
COLWI CK  array . 

Figure  13.  High  resolution  vertical  compoment  acceleration  wavenumber 
power  spectra  at  3*7  Hz  for  the  first  5.12  sec  of  data. 

The  origin  corresponds  to  the  slowness  of  the  first  arrival. 

The  azimuth  of  the  source  is  in  the  direction  of  the  negative 

Y  axis. 

Figure  l^*.  Similar  to  Figure  13  for  the  transverse  component. 

Figure  15.  Example  of  a  cross-correlation  after  normalization  by  the 
corresponding  auto-correlat ions. 

Figure  16.  Contours  of  the  spatial  dependence  of  the  maximum  of  the 
two-station  cross-correlations.  The  dots  represent  the 
postions  of  the  13  independent  data  points  which  constrain 
the  surface.  Values  at  the  data  points  are  suppressed 


for  clarity. 
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Figure  17- 


Figure  l8. 


Figure  19- 


Spatial  decay  of  the  two-station  cross-correlations.  The 
curves  represent  sections  through  the  surface  of  Figure  16: 
solid  -  direction  of  maximum  gradient,  short  dashed  -  transverse 
direction,  long  short  dashed  -  longitudinal  direction. 

Example  of  a  normalized  and  low-pass  filtered  cross-correlation 
function  at  corner  frequencies  of  2.5,  5-0,  10.0,  and  20.0  Hz. 
Filtered  cross-correlations  versus  corner  frequency.  Seven  pairs 
of  stations  are  presented,  either  pure  transverse  separation 
or  pure  longitudinal  separation. 
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COLWICK  ARRAY  WAVENUMBER  IMPULSE  RESPONSE 
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FILTERED  TWO  STATION  CROSS  -  CORRELATIONS 
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III.  REGIONAL  CRUSTAL  MODEL  FROM  SURFACE  WAVES 

In  previous  reports,  we  discussed  the  proposed  seismic  studies  using 
seismograms  recorded  at  regional  distances  (0  -  10°),  summarizing  the  type 
of  surface  wave  data  available  for  analysis,  the  planned  method  of  computing 
observed  dispersion  and  attenuation  data,  and  a  specific  planned  crustal 
investigation  at  the  Nevada  Test  Site  (NTS). 

Through  a  substantial  effort  in  the  past  year,  we  have  succeeded  in  the 
development  of  methods  to  analyze  and  interpret  a  large  body  of  regional 
broadband  surface  wave  data  generated  by  nuclear  explosions  and  collapses 
at  NTS.  Following  Herrin  and  Goforth  (1977),  we  developed  a  Phase-Matched 
Filter  (PMF)  Program  and  other  computer  routines  to  be  used  in  the  problem 
of  obtaining  Q  and  velocity  structure  from  observed  group  velocity,  phase 
velocity,  and  attenuation  coefficients  of  fundamental  and  higher  mode  surface 
waves. 

Numerical  tests  show  that  PMF  is  superior  to  filtering  techniques  such 
as  Moving  Window  Analysis  (MWA)  of  Landisman,  Dziewonski  and  Sato  (1969). 

In  M\'/A,  modal  interference  and  large  secondary  arrivals  often  cause  observed 
dispersion  curves  to  be  discontinuous.  In  PMF,  however,  we  take  advantage 
of  the  approximately  known  regional  dispersion  characteristics  to  construct 
a  filter  which  we  apply  to  digitized  surface  wave  data.  As  a  result,  the 
signal  dispersion  is  removed.  In  the  time  domain,  it  corresponds  to  the  com¬ 
pression  of  signal  energy  to  a  short  symmetrical  pulse.  Random  noise  and 
multipathing  energies  are  removed  by  appropriate  windowing  centered  at  the 
pulse  maximum.  Iterative  application  to  seismograms  recorded  at  several 
sites  yields  both  group  velocity  and  attenuation  properties  of  the  uncon¬ 


taminated  waveform. 
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Figures  1  and  2  show  a  typical  result.  Shown  are  the  seismogram  and 
two  sets  of  observed  Rayleigh  wave  dispersion  data  obtained  by  applying  PMF 
and  MWA,  respectively,  to  the  broadband  data  recorded  by  the  seismic  network 
of  the  Lawrence  Livermore  National  Laboratory.  The  filter  used  in  the  PMF 
has  been  constructed  from  a  modified  velocity  model  of  the  Great  Basin  by 
Priestly  and  Brune  (1978) .  We  found  that  MWA  failed  to  yield  useful  dis¬ 
persion  data  outside  the  period  range  7-30  seconds.  Application  of  PMF  to 
digitized  surface  wave  data  generated  by  large  events  (mag  ^  5.5)  in  the  dis¬ 
tance  range  of  2-4“  usually  produce  excellent  dispersion  and  amplitude  data 
in  the  period  range  5-50  seconds,  which  represent  a  substantial  broadening 
of  the  data  bandwidth.  Confidence  in  the  extended  data  came  from  the  fact 
that  results  for  a  given  station  are  consistent  from  event  to  event,  but 
different  from  values  at  other  stations  and  from  the  starting  model  (represented 
by  the  input  group  velocity  curve) . 

Figures  3  and  4  show  the  partial  derivatives  of  Rayleigh  wave  phase 
velocity  and  the  surface  wave  Qr  with  respect  to  shear  velocity  for  the 
slightly  dissiptative  continental  crustal  model  of  Table  1.  The  result 
illustrate  that  amplitude  data  can  play  a  significant  role  in  surface  wave 
inversion  because  they  are  strongly  dependent  on  the  velocity  structure,  and 
also  that  the  high  sensitivity  of  the  phase  velocity  and  amplitude  of  the 
first  higher  mode  to  the  crustal  shear  wave  velocity  offers  a  means  to  Improve 
depth  resolution.  Emphasis  is  now  on  experimental  tests  to  find  an  approximate 
strategy  for  inversion  and  starting  model  selection  for  the  NTS  data  set. 
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TABLE  1 


Depth 

km 

V 

P 

km/ sec 

V 

s 

km/ sec 

% 

% 

Density 

gm/cc 

0.0 

5.20 

3.30 

600 

300 

2.60 

6.0 

6.40 

3.90 

1000 

500 

3.00 

36.0 

8.00 

4.62 

1600 

800 

3.35 

Crustal  model  used  for  data  shox^m  in  Figures  3  and  4. 
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IV.  NEAR-SOURCE  EFFECTS  ON  P  WAVES 
(A  review  submitted  in  response  to  a  DARPA  request  in  February  I98O) 

The  matter  of  the  generation  of  elastic  waves  by  an  explosion 
buried  In  the  earth  would  appear  to  be  a  relatively  simple  problem. 

Thus  the  exercise  of  characterizing  an  explosion  through  the  analysis 
of  the  radiated  elastic  waves  should  be  quite  tractable  and  yield 
definitive  results.  In  practice,  this  has  not  proved  to  be  a  simple 
process,  and  there  remain  several  unanswered  questions  concerning  the 
relationship  between  an  explosive  source  and  the  waves  that  are  emitted 
from  the  source  region. 

Part  of  the  difficulty  undoubtedly  relates  to  the  fact  that  observed 
seismic  waves  contain  the  combined  effects  of  the  source  and  the 
propagation  between  source  and  receiver,  and  it  is  often  difficult  to 
separate  these  two  effects.  However,  with  current  knowledge  about 
earth  structure  and  the  improving  capability  for  calculating  the  propagation 
effects,  only  a  diminishing  amount  of  the  difficulty  can  be  attributed 
to  this  cause.  In  what  follows,  therefore,  it  will  be  assumed  that 
propagation  effects  outside  the  immediate  source  region  can  be  calculated 
and  are  not  a  major  obstacle  in  the  interpretation.  However,  it  is 
important  to  keep  in  mind  this  basic  fact  that  the  effects  of  source  and 
propagation  are  often  indistinguishable,  and  unless  one  is  known  the 
other  can  not  be  uniquely  determined. 

Setting  aside  the  effects  of  propagation  outside  the  source  region, 
we  must  consider  the  possibi 1 ities  that  the. waves  generated  by  an 
explosion  are  not  as  simple  as  we  might  expect  or  that  the  waves  are 
modified  by  effects  very  near  the  source.  In  most  cases  it  is  not  very 
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meaningful  to  try  to  separate  these  two  possibilities,  and  so  it  is 
useful  to  combine  them  and  consider  the  Waves  that  propagate  outward 
from  a  general  source  region,  which  is  taken  to  be  a  region  surrounding  the 
explosion  and  including  part  of  the  crust  in  the  immediate  vicinity. 

The  subject  of  this  summary  will  be  our  current  understanding  of 
these  waves  that  propagate  outward  from  the  source  region  of  a  buried 
explosion.  The  discussion  will  be  mainly  restricted  to  the  P  waves, 
and  thus  only  body  waves  will  be  considered. 

Theoretical  Considerations.  The  problem  which  has  the  most  physical 
similarity  to  that  of  a  buried  explosion  and  also  has  an  analytic  closed-form 
solution  is  that  of  a  pressure  pulse  applied  to  the  interior  of  a  spherical 
cavity  in  a  homogeneous  elastic  medium.  Although  many  differences  remain 
between  this  idealized  mathematical  problem  and  the  actual  situation  of 
an  explosion  in  the  earth,  it  is  reasonable  to  expect  that  its  solution 
might  provide  a  first  approximation  to  the  more  complicated  problem. 

The  solution  to  the  problem  of  a  pressure  pulse  in  a  spherical 
cavity  is  well-known,  and  numerous  treatements  can  be  found  in  the 
literature  (for  example:  Jeffreys,  1931;  Sharpe,  19^2;  Blake,  1952; 

Favreau,  1969)-  Because  of  the  spherical  symmetry  of  the  problem,  the 
solution  can  be  expressed  in  terms  of  a  scalar  function  known  as  the 
reduced  displacement  potential,  which  is  only  a  function  of  the  reduced 
travel  time.  The  usual  procedure  is  to  specify  the  pressure  history 
within  the  cavity,  solve  for  the  reduced  displacement  potential,  and 
then  obtain  the  displacement  at  any  point  by  taking  the  gradient  of 


the  reduced  displacement  potential  divided  by  the  distance  from  the 
source. 

Because  of  its  simple  form,  the  reduced  displacement  potential 
has  become  a  popular  device  for  characterizing  an  explosive  source. 

Since  it  is  independent  of  distance,  it  can  be  determined  at  any 
convenient  distance  from  the  source.  This  is  helpful  in  dealing  with 
the  complications  of  an  inelatic  zone  which  surrounds  large  explosions. 
The  initial  cavity  of  a  contained  explosion  will  in  general  be  surrounded 
by  successive  zones  of  vaporization,  melting,  cracking,  and  inelastic 
stresses  before  reaching  a  zone  where  the  assumptions  of  linear  elasticity 
are  appropriate.  In  place  of  the  radius  of  the  initial  cavity,  it  is 
customary  to  use  the  inner  boundary  of  the  region  of  elastic  behavior 
as  the  effective  source  radius.  This  is  often  called  the  elastic  radius. 
With  this  modification  of  the  problem  the  pressure  history  within  the 
cavity  gets  replaced  by  the  stress  wave  which  arrives  at  the  elastic 
radius. 

It  should  be  noted  in  passing  that  the  problem  where  the  source 
consists  of  shear  stresses  applied  to  the  interior  of  a  spherical 
cavity  also  has  known  solutipns  (Jeffreys,  1931;  Honda,  I960).  However, 
so  far  these  results  seem  to  have  found  little  application  in  the  problem 
of  a  buried  explosion. 

A  spherically  symmetric  source,  such  as  a  pressure  pulse  in  a 
spherical  cavity,  generates  only  P  waves.  However,  if  such  a  source  is 
placed  in  a  medium  with  a  preexisting  shear  stress,  then  S  waves  will 
also  be  generated  (Archambeau,  1972).  Closely  related  to  this  mechanism 
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is  that  of  triggering  a  tectonic  earthquake  on  a  nearby  fault  (Andrews, 
1973).  In  either  case  the  secondary  source  related  to  the  shear  stress 
has  the  form  of  a  shear  dislocation,  and  so  its  radiation  pattern  is  a 
quadrapoie  and  It  is  roughly  an  order  of  magnitude  more  efficient  in 
generating  S  waves  than  P  waves.  Thus,  observational  studies  of  this 
effect  have  been  based  primarily  on  S  waves  and  surface  waves  (Press 
and  Archambeau,  1962;  Aki  et  al.,  1969;  Archambeau  and  Sammis,  1970; 

Toksoz  efal.,  1971;  Lambert  et  al.,  1972;  Aki  and  Tsai,  1972;  Toksoz 
and  Kehrer,  1972).  It  appears  from  the  results  of  these  studies  that, 
while  reasonable  levels  of  prestress  can  be  Important  in  the  generation 
of  S  waves  and  surface  waves,  the  direct  P  waves  generated  by  this  secondary 
source  will  be  small  compared  to  those  generated  by  the  explosive  source. 
However,  the  mechanism  whereby  S  waves  from  the  secondary  source  are 
converted  to  P  waves  at  nearby  boundaries  could  be  an  important  source 
of  P  waves. 

The  next  step  in  constructing  a  more  realistic  model  for  a  buried 
explosion  in  the  earth  is  to  consider  the  effects  of  material  inhomogeneity 
in  the  vicinity  of  the  explosion.  Such  inhomogeneity  causes  the  reflection 
and  refraction  of  primary  waves  from  the  source,  the  conversion  of  wave 
type  from  P  to  S  and  vice  versa,  and  the  generation  of  interface  waves. 

For  shallow  explosions,  a  very  important  inhomogeneity  is  the  free  surface 
of  the  earth.  Geologic  layering,  the  water  table,  and  fault  surfaces 
are  other  types  of  inhomogeneities  which  can  also  have  significant  effects. 

Attempts  to  obtain  exact  solutions  to  the  problem  of  an  explosion 
in  the  vicinity  of  an  inhomogeneity  have  not  been  very  successful. 


-  43  - 


The  problem  of  a  finite  spherical  source  embedded  in  a  homogeneous 
elastic  halfspace  has  been  considered  by  Ben-Menahem  and  Cisternas  (1963) 
and  by  Thi ruvenkatachar  and  Viswanathan  (1965,  1967).  This  is  a  problem 
with  mixed  boundaries  and  its  solution  is  very  difficult,  the  answer 
usually  being  expressed  as  an  infinite  series.  Because  of  the  rather 
untractable  form  of  the  results,  the  analytical  treatment  of  this  problem 
has  not  yet  contributed  any  practical  results  to  the  problem  of  a 
buried  explosion. 

Provided  one  is  willing  to  approximate  an  explosion  with  a  point 
source,  the  effect  of  vertical  inhoraogenei ty  in  the  vicinity  of  the 
source  can  be  handled  in  a  satisfactory  manner.  Two  approaches  are 
commmonly  used,  the  method  of  generalized  rays  (Helmberger  and  Harkrider, 
1972)  or  the  method  using  propagator  matrices  (Fuchs,  1966;  Hasegawa, 
1971). 

The  treatment  of  lateral  inhomogeneity  near  an  explosive  source 
is  a  more  difficult  problem.  This  would  include  such  features  as  dipping 
layers,  faults,  and  surface  topography.  In  general,  numerical  methods 
are  required  to  calculate  the  effects  of  this  type  of  inhomogeneity, 
but  usually  the  detailed  knowledge  about  the  geometry  of  such  features 
is  not  sufficient  to  justify  an  elaborate  computational  treatment. 

The  presence  of  inhomogeneity  near  the  explosive  source  can  also 
lead  to  additional  inelastic  effects.  Spallation  is  an  important  effect 
of  this  type.  When  the  P  wave  from  an  explosion  is  reflected  at  the 
free  surface  with  a  change  in  sign,  the  associated  stress  can  cause 
failure  in  tension  of  the  near  surface  material.  When  failure  occurs. 
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part  of  the  near-surface  material  may  separate  and  move  bal 1 i sti cal 1y 
upward,  eventually  falling  back  and  impacting  the  earth  at  some  later 
time.  This  closure  of  the  spall  is  sometimes  referred  to  as  slapdown. 

This  process  of  spallation  has  been  fully  described  and  documented  in 
the  literature  (Eisler  and  Chilton,  1964;  Eisler  et  al.,  1966;  Chilton 
et  al.,  1966;  Perret,  1972;  Viecelli,  1973;  Springer,  1974).  The  energy 
in  the  P  wave  which  leaves  the  source  in  an  upward  direction  is  converted 
into  a  reflected  pP  wave,  a  reflected  sP  wave,  a  surface  wave,  inelastic 
effects,  and  the  waves  generated  by  the  slapdown.  Viewed  from  a  distance, 
the  primary  effects  of  spallation  upon  the  P-wave  coda  is  a  diminished 
amplitude  of  the  pP  wave  and  an  additional  phase  at  the  time  of  slapdown. 

On  the  basis  of  the  preceding  discussion,  it  is  possible  to  construct 
a  general  model  of  a  buried  explosion.  The  mode)  is  primarily  an  elastic 
model  so  it  begins  with  a  stress  pulse  applied  at  the  elastic  radius. 

This  generates  an  outward  propagating  P  wave  and,  if  the  source  region  is 
prestressed,  also  an  S  wave.  These  outward  propagating  waves  interact 
with  inhomogenei t ies  in  material  properties  near  the  source  to  produce 
reflected  and  interface  waves.  Waves  reflected  from  above  the  source  may 
also  cause  spallation  and  the  associated  slapdown.  The  outward  propagating 
waves  may  also  trigger  secondary  shear  dislocations  at  near-by  stress 
concentrations.  All  of  these  effects  combine  and  interact  to  produce  the 
waves  that  propagate  out  from  the  general  source  region. 
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Experimental  Results.  Consider  now  the  question  of  what  can  be 
learned  about  a  buried  explosion  from  an  analysis  of  experimental  data. 

This  is  a  basic  inverse  problem.  The  effects,  primarily  the  waveforms 
of  elastic  waves,  are  to  be  used  to  estimate  properties  of  the  cause, 
the  explosive  source.  The  usual  approach  is  to  construct  a  general 
model  of  the  source  containing  a  number  of  undetermined  parameters. 

The  observational  data  are  then  analyzed  to  deterimine  whether  the 
model  is  capable  of  explaining  the  data  and,  if  so,  what  values  should 
be  given  to  the  parameters. 

Because  of  the  advantages  already  mentioned,  it  has  become  common 
to  characterize  an  explosive  source  by  Its  reduced  displacement  potential, 
and  so  the  experimental  determination  of  the  reduced  displacement  potential 
has  received  considerable  attention.  The  generalized  form  of  the  reduced 
displacement  potential  as  a  function  of  time  is  an  abrupt  start,  a  smooth 
rise  to  a  maximum  over  a  finite  time,  and  then  a  decrease  to  a  static 
level.  Three  numbers  -  the  rise  time,  the  static  value,  and  the  ratio  of 
maximum  to  static  values  -  describe  the  major  features  of  such  a  function. 

In  the  frequency  domain,  the  first  time  derivative  of  the  reduced  displacement 
potential  has  a  correspondingly  simple  form.  Going  from  high  to  low 
frequencies,  the  spectrum  rises  at  some  slope,  reaches  a  maximum  near  what 
is  called  the  corner  frequency,  and  then  decreases  to  a  constant  value 
at  low  frequencies.  In  terms  of  the  time  domain  parameters,  the  spectral 
amplitude  scales  with  the  static  value,  the  frequency  scales  with  the 
rise  time,  the  ratio  of  the  spectral  maximum  to  the  low  frequency  level 
depends  upon  the  ratio  of  maximum  to  static  values  of  the  reduced  displacement 
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potential,  and  the  high  frequency  slope  depends  upon  the  abruptness  of 
the  beginning  of  the  reduced  displacement  potential. 

The  simple  model  of  a  pressure  pulse  applied  to  the  interior  of  a 
spherical  cavity  can  be  used  to  relate  the  parameters  of  the  reduced 
displacement  potential  to  the  physical  properties  of  the  explosive 
source.  The  static  value  depends  upon  the  static  pressure,  the  cavity 
volume,  and  the  material  properties.  The  rise  time  depends  upon  the 
cavity  radius  and  the  material  properties.  The  ratio  of  maximum  to 
static  values,  sometimes  called  the  overshoot,  depends  upon  the  time 
history  of  the  pressure  pulse  and  the  material  properties.  Considering 
the  pressure  pulse  to  consist  of  two  main  parts,  an  impulse  and  a  step, 
the  overshoot  will  increase  as  the  ratio  of  impulse  to  step  increases. 

For  a  more  complete  parameterization  of  the  reduced  displacement 
potential,  it  can  be  approximated  with  an  analytic  function.  Haskell  (196?) 
argued  that  displacement,  velocity,  and  acceleration  should  all  be 
continuous  at  the  elastic  radius  and  used  a  fourth  order  polynomial  in 
time.  Von  Seggern  and  Blandford  (1972)  required  that  only  the  displacement 
be  continuous  and  used  a  second  order  polynomial.  Mueller  and  Murphy 
(1971)  used  the  theoretical  solution  for  a  pressure  pulse  within  a  cavity 
and  a  semi-empirical  expression  for  the  shape  of  the  pressure  pulse  to 
arrive  at  an  analytic  expression  for  the  reduced  displacement  potential. 

The  basic  difference  in  these  three  models  is  primarily  at  the  high 
frequencies,  where  the  Haskell  model  falls  off  with  a  -4  slope  on  a  log-log 
scale,  while  both  the  von  Seggern-Blandford  and  Muel ler-Murphy  models 
fall  off  with  a  -2  slope.  The  experimental  data  (von  Seggern  and  Blandford, 
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1972;  Murphy,  1977;  Burdick  and  Helmberger,  1979)  at  teleseismic,  regional, 
and  near  distances  mostly  favor  the  -2  slope  of  the  von  Seggern-Blandford 
and  Muel ler-Murphy  models.  From  an  analysis  of  very-near  data,  Peppin  (1976) 
found  evidence  for  a  slope  of  at  least  -3  at  the  high  frequencies. 

The  problem  of  estimating  the  reduced  displacement  potential  for  an 
explosion  can  be  approached  from  at  least  three  directions.  The  dynamic 
equations  for  the  explosion  and  the  surrounding  inelastic  region  can  be 
solved  numerically  and  the  calculations  carried  out  to  the  elastic 
radius  (Holzer,  1966;  Rodean,  1971).  Another  approach  is  to  measure 
the  ground  motion  as  near  tlie  source  as  possible  while  still  in  the 
elastic  region,  and  then  calculate  the  reduced  displacement  potential 
on  the  basis  of  these  measurements  (Werth  et  al.,  1962;  Werth  and 
Herbst,  1963).  In  most  cases  the  reduced  displacement  potential  can 
only  be  determined  for  times  less  than  about  0.5  sec  because  the  effect 
of  the  free  surface  and  other  departures  from  spherical  symmetry  begin 
to  affect  the  results  at  later  times.  A  third  approach  is  to  record 
elastic  waves  at  near  to  teleseismic  distances  and  then  attempt  to 
infer  the  reduced  displacement  potential  which  best  explains  the 
observations.  What  emerges  from  this  approach  is  an  apparent  reduced 
displacement  potential,  because,  as  already  discussed,  the  waves  that 
emerge  from  the  general  source  region  can  consist  of  considerably  more 
than  the  direct  P  wave  from  the  explosion.  Such  additional  effects 
must  be  taken  into  account  in  the  interpretation  of  the  data. 

Because  of  the  abundance  of  easily  accessible  data,  numerous  studies 


have  used  waveform  data  recorded  at  near  to  teleseismic  distances  to 


investigate  the  details  of  explosive  sources.  Various  methods  of 
interpretation  have  been  employed.  A  few  of  the  representative 
studies  will  be  summarized  below. 

Molnar  (1971),  Kulhanek  (1971),  and  Wyss  et  al.  (1971)  all 
studied  the  spectra  of  teleseismic  P  waves  and  found  that  the  spectra 
were  modulated.  This  can  be  explained,  at  least  partly,  by  the 
interference  of  the  P  and  pP  waves.  Filson  and  Frasier  (1972)  and 
King  et  al.-  (1972)  fit  theoretical  models  to  the  spectra  of  teleseismic 
P  waves  to  estimate  parameters  of  the  reduced  displacement  potential. 

Aki  et  al.  (197^)  studied  both  near  and  teleseismic  data  and 
found  evidence  for  a  large  overshoot  ratio  in  the  reduced  displacement 
potential.  Peppin  (1976)  analyzed  near  and  regional  data  and  did  not 
find  evidence  for  overshoot.  Murphy  (1977)  studied  a  variety  of  data 
from  near  to  teleseismic  distances  and  concluded  that  the  data  were 
consistent  with  the  source  model  of  Mueller  and  Murphy  (1971). 

Frasier  (1972)  deconvolved  teleseismic  P  waves  and  found  evidence 
of  a  pP  phase  plus  a  later  phase,  possibly  due  to  slapdown.  Bakun 
and  Johnson  (1973)  applied  homomorphic  deconvolution  to  teleseismic 
P  waves  and  found  indications  of  the  pP  reflection  and  also  a  later 
phase  which  was  consistent  with  an  interpretation  in  terms  of  slapdown. 
Burdick  and  Helmberger  (1979)  used  synthetic  seismograms  to  model 
teleseismic  P  waves  and  found  that  the  data  could  be  explained  by 
substantial  overshoot  in  the  reduced  displacement  potential  and 
reflected  crustal  phases,  but  did  not  require  a  slapdown  phase. 
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Stump  and  Johnson  (1977)  and  Stump  (1979)  have  developed  a  general 
inverse  method  for  estimating  the  second-rank  seismic  moment  tensor. 

The  trace  of  this  tensor  is  equivalent  to  the  reduced  displacement 
potential  and  the  deviatoric  components  provide  a  means  of  expressing 
other  effects,  such  as  tectonic  stress  release. 

There  is  obviously  not  total  agreement  among  the  observational 
studies  concerning  the  source  properties  of  explosions.  This  is  partly 
due  to  the  different  methods  of  interpretation  which  have  been  employed. 

It  is  also  due  to  a  certain  degree  of  nonuniqueness  that  exists  in  the 
basic  problem.  This  is  compounded  by  the  fact  that  propagation  effects 
must  usually  be  taken  into  account  in  the  interpretat ion  of  the  data, 
and  any  uncertainty  in  earth  structure  can  get  translated  into  nonuniqueness 
in  the  source  properties.  For  instance,  a  peaked  spectrum  of  the  P 
wave  coda  can  be  produced  in  at  least  three  ways:  overshoot  of  the 
reduced  displacement  potential,  interference  caused  by  reflected  waves, 
or  ovei — correction  for  the  effects  of  attenuation  in  the  interpretation 


process. 
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Cone  1  us i ons .  The  available  theoretical  models  of  an  explosive  source 
embedded  in  a  realistic  crustal  structure  appear  to  be  sufficiently  general 
to  explain  the  major  features  of  the  observational  data.  It  is  clear  that 
considerably  more  than  the  direct  P  wave  emerges  from  the  general  source 
region.  Reflections  from  inhomogeneities  near  the  explosion,  particularly 
the  free  surface,  are  definitely  important.  i.ie  importance  of  spallation 
and  slapdown  upon  the  generation  of  P  waves  is  still  rather  uncertain. 
Tectonic  strain  release  may  be  important  in  some  instances,  and,  when 
it  is  significant,  it  is  probably  not  the  direct  P  waves  from  the  secondary 
sources  but  rather  the  S  to  P  converted  phases  which  contribute  most  to 
the  P  wave  coda.  The  entire  matter  of  how  S  waves  are  generated  in  the 
vicinity  of  an  explosion  is  still  not  completely  understood,  and  some 
new  process,  such  as  the  acoustic  fluidization  suggested  by  Nelosh  (1979), 
may  eventually  provide  the  answer. 

The  methods  of  interpreting  the  observational  data  to  infer  source 
properties  are  steadily  improving,  and  the  increased  use  of  synthetic 
seismograms  and  more  complete  inversion  schemes  should  be  very  helpful. 

This  prog''ess  in  interpretation  will  require  a  more  accurate  knowledge 
of  earth  sturcture,  because  it  is  doubtful  that  the  accuracy  of  inferred 
source  prope'^ties  can  ever  be  greater  than  that  of  the  earth  model  used 
in  the  interpretation.  The  anelastic  properties  of  the  earth  present 
a  current  problem  in  this  respect,  because  the  correction  for  attenuation 
which  must  usually  be  applied  in  the  int*".  pretation  process  is  still 
rather  uncertain  and  this  can  have  a  major  effect  at  the  high  frequencies. 
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Hopefully,  continued  improvements  in  the  methods  of  interpretation  and 
knowledge  of  earth  structure  will  reduce  the  inherent  nonuniqueness 
in  this  inverse  problem. 

The  reduced  displacement  potential  has  become  a  popular  means  of 
summarizing  the  properties  of  an  explosive  source  and  has  been  quite 
useful.  But  is  should  be  emphasized  that  in  most  cases  this  is  only 
an  apparent  reduced  displacement  potential.  The  requirement  for  spherical 
symmetry  holds,  if  at  all,  only  for  a  few  tenths  of  a  second,  and 
asymmetry  in  the  source  region  is  an  important  factor  in  most  P  wave 
codas.  The  importance  of  this  point  is  the  realization  that  the 
apparent  reduced  displacement  potential  for  a  given  event  may  be 
different  as  viewed  from  different  distances,  different  azimuths,  or 
different  types  of  data. 
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M .  Archival  of  Digital  Seismic  Data 

A  catalog  has  been  assembled  for  most  of  the  digital  data  which  we 
have  collected  over  the  years  in  the  process  of  conducting  research  sponsored 
by  DARPA.  Each  entry  in  the  catalog  consists  of  a  set  of  data  which  have 
been  obtained  by  digitizing  a  portion  of  a  seismogram  or  a  suite  of  seismograms 
that  were  recorded  from  a  particular  seismic  event.  Most  of  the  events  in 
the  catalog  are  explosions  or  earthquakes  recorded  at  local  or  regional 
distances  by  seismograph ic  stations  in  California,  Nevada,  and  Utah,  although 
a  few  events  recorded  at  teleseismic  distances  are  also  included.  The 
present  catalog  contains  about  200  entries  arranged  in  chronological  order 
for  the  time  period  August  I965  to  November  1971. 

The  digital  data  are  now  stored  on  punched  cards  in  a  binary  format, 
and  the  possibility  of  putting  the  data  on  digital  magnetic  tape  is  being 
considered.  Plans  call  for  the  time  period  covered  by  the  catalog  to  be 
extended  from  1971  up  to  the  present.  We  also  plan  to  put  the  index  to 
the  catalog  in  a  computer  where  it  can  be  accessed  by  a  text  editor  and 
thus  ease  the  task  of  updating  the  index. 

The  following  pages  contain  the  first  version  of  the  catalog  index. 

It  contains  the  information  that  was  readily  available  in  a  first  pass 
through  the  raw  data  files.  Considerable  more  information  about  the  digital 
data  is  available  in  various  reports  and  records. 


-  57  - 


What  follows  is  a  brief  explanation  of  the  information  which  is 
contained  in  the  index  of  the  digital  data  catalog. 

DATE  -  Date  of  the  seismic  event. 

EVENT  -  Type  of  seismic  event. 

EQ  -  Earthquake. 

EX  -  Nuclear  explosion. 

C  -  Collapse  of  a  nuclear  explosion. 

EXA  -  Event  following  a  nuclear  explosion  which  is  not 
identified  as  a  collapse. 

BLAST  -  Chemical  explosion. 

STA.  -  Seismic  station  which  recorded  the  event. 

JAS  -  Jamestown,  California  (U.C.  Berkeley  network). 

BRK  -  Berkeley,  California  (U.C.  Berkelely  network). 

PRI  -  Priest,  Cal ifornia  (U.C.  Berkeley  network). 

MMG  -  Mina,  Nevada  (LLNL  network). 

EMG  -  Elko,  Nevada  (LLNL  network). 

KMG  -  Kanab,  Utah  (LLNL  network). 

LMG  -  Landers,  California  (LLNL  network). 

BEFL  -  Belleview,  Florida  (LRSM  network). 

HNME  -  Moulton,  Maine  (LRSM  network). 

Kanab,  Utah  (LRSM  network). 
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LCNM  -  Las  Cruces,  New  Mexico  (LRSM  network). 

RKON  -  Red  Lake,  Ontario  (LRSM  network). 

SJTX  -  San  Jose,  Texas  (LRSM  network). 

TIME  -  Approximate  origin  time  of  event  (hours,  minutes  GMT). 

INST.  -  Type  of  instrument  used  to  record  the  event. 

BBV  -  Broad  band  velocity  (see  Peppin  and  McEvilly,  Geophys. 

J.  R.  Astr.  Soc.,  37,  227-243,  1974). 

SP  -  Short  period. 

FBA  -  Force-balance  accelerometer. 

COMPONENTS  -  Component  of  ground  motion  and  other  information. 

Z  -  Vertical. 

L  -  Longitudinal  with  respect  to  source  direction. 

T  -  Transverse  with  respect  to  source  direction. 

SM  -  Strong  motion,  low  gain. 

P  -  P  wave  arrival . 

Z  -  Leader  (recorded  immediately  preceding  the  event). 

COMMENTS  -  Any  other  available  information  about  the  event. 
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